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1. Introduction 


There have been a number of recent papers estimating high-energy 
neutrino fluxes and spectra from various astrophysical processes (Berezinsky 
and Zatsepin 1970, Stecker 1973a, Berezinsky and Zatsepin 1977 and references 
therein, Silberberg and Shapiro 1977, Margolis, Schramm and Silberberg 1978, 

Eichler 1978.) In this work, I reexamine the problem by 1) presenting 

the results of a detailed calculation of galactic v~production in cosmic-ray 

interactions with interstellar gas and extragalactic v-production by inter- 
actions of ultrahigh energy cosmic rays with the 3K universal background 
radiation, 2) using these results together with present y“ ra y observations 
to predict neutrino fluxes and event rates, and 3) reexamining some of the 
assumptions and results of previous work, pointing out significant differences, 
2, Production Rates 

The first basic production process for high energy cosmic neutrinos 
is the deca} of charged pions produced in cosmic ray interactions with 
interstellar gas, primarily pp interactions. This process will henceforth 
be referred to as "pp", although the effects of ap pHe and aHe interactions 
are included in the calculation. The second process involves the photo- 
production of TT-mesons by interaction of ultrahigh energy cosmic rays with 
the 3K universal microwave background radiation (henceforth referred to as 
yp) and subsequent meson decay. Both of these types of interactions involve 
the accompanying production of tt° mesons and their decay Into cosmic y-rays. 
Thus, the production rates of cosmic y-rays and neutrinos are generically 
linked. 

A detailed discussion of the kinematics of the production and decay of 
sec^dary particles produced in "pp" interactions may be found in Stecker 
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(1971) • Details of the yp process have also been previously given (Stecker 
1968). The ^-ray production spectrum In yp Interactions has been calculated 
by Wdowczyk et al, (1972), Stecker and Bredekamp (1972), and Stecker (1973a). 

In the present paper, the "pp" itrino production spectrum was 

calculated for pp Interactions up to 30 GeV by methods previously employed 
(Stecker 1971, 1973b). At higher energies, scaling was assumed to hold. At these 
energies, the parameters of Ganguli and Sreekantan (1976) were adopted for 
the rapidity distribution of charged plons. These authors have calculated 
y-ray spectra which are In good agreement with those of Stecker (1970) for 
y-ray energies < 10 GeV, however, at higher energies, the assumption of 
scaling gives larger fluxes and a flatter spectrum than those calculated by 
stacker (1971) using an "Isobar + fireball" model in "pre-scallng days". 

It should be noted that the isobar + fireball (I-F) model is equivalent 
to a "leading-pion" model at high energies, since the Isobar carries off 
~ 50% of the energy and decays into a "leading plon". Both scaling and 
the I-F models produce a secondary spectrum which has the same spectral 
index as the primary spectrum in the high-energy limit. (See Appendix I). 

The results of the ' pp" production spectrum calculation are shown in 
Fig. 1. The two neutrino production spectra are given for an Interstellar 

. p 

hydrogen density of 1cm ~ so that it is really a production rate per hydrogen 
atom. The upper neutrii o curve and the y-ray production spectrum shown 


In the figure are calcu ated using a primary cosmic ray spectrum 

t (E) = 2.35 E" 2o67 cm" 2 r 1 sr“ 1 GeV” 1 forEp>10GeV and the lower curve is 

“2 75 "2 -1 —1 

vbtalued for a primary -pectrum I^CEp ) * 2,0 E 0 cm s sr GeV (Ryan, 


;*■ me aid Balasubrahraanyan 1972), The y-ray spectrum is for all y-rays 
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from n*-decay; the v-spectrum Is for the v„ component. Since each plon decay 

result* in one one v^and one v e (v e )» all of roughly the same energy. 

the tdlal number of v's produced (although of four different types) Is a 

factor of 3 higher than that shown for v^’s. In this regard, it should be 

kept in mind that the cross section ratios given in Table 1 apply. 

The production rate for yP interactions has been calculated using 

the method of Stecker and Bredekamp (1972) o As an example of the range 

of uncertainty, two different primary spectra were used, i.e„ I p (Ep) ■ 

6.4xl0~ 31 (E p /10 11 )" 2,6 cm“ 2 s" 1 sr" 1 GeV" 1 (Linsley 1963) and Ip(E p ) - 

2.4xlO~ 31 (E /10 11 )“ 3o24 cm“ 2 s~ 1 sr” 1 GeV“ 1 (Andrews et al. 1971). The 
P 

latest analysis of the ultrahlgh energy air-shower data (Hillas, private 

communication) gives Ip(E p ) «*• 4x10 3 ^(E p /10^) 3 cm 2 s ^sr ^GeV in the 
8 10 

energy range 10 ^E p <L0 GeV with indications of a flattening to I p (E p ) 

_-»n 11 .9 i in 11 

2x10 (E p /10 A ) 0 in the energy range 10 <E p <L0 GeV. The results of 

the yP calculation are shown in Figure 2 for the two spectra chosen. The 
right-hand scale of Figure 2 also shows the diffuse background flux from 
this process obtained by multiplying by the factor c/(4ttH q ) where H q , the 
Hubble constant, is taken to be 50km s Hlpc \ This flux, whi"h is alro 
shown in the subsequent figures, only holds if the ultrahigh energy primary 
spectrum is universal, an assumption which is contradicted by the lack of 
an observed high energy cutoff in the spectrum (Greisen 1966, Zatsepin and 
Kuzmin 1966, Stecker 1968, 1978a). Thus the fluxes given may be over- 
estimated and may actually be upper limits. 

Figure 3 shows the integral galactic ’’pp" neutrino production spectrum. 
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3. Diffuse Cosmic Neutrino Fir 

Figures 4 and 5 show the dinerential and integral neutrino fluxes from 
various sources. In figure 4, the croashatched region marked yp is obtained from 
the curves shown in Figure 2. The hatched region marked pp(G,C, |b |*jlO°) is 
for galactic v-production coming from the galactic central region defined 
by galactic longitude 330°s-t^40° and latitudes Jb|^10° o The dashed line is 
the flux computed using the isobar-fireball (I-F) model (Stecker 1971), which 
is probably too low. The lines bounding the hatched region are obtained 
using the two ”pp" production spectra shown in Figure 1, All galactic 
"pp" curves are noramlized by assuming that 1/2 of the galactic y-radiation 
above 100 MeV observed by SAS-2 (Fichtel et al 0 1978) is from TT°-decay 
and by relating the y-ray and v production using the results given in 
Figure 1„ Spectral measurements indicate that ~ 50% of the galactic y-rays 
observed above 100 MeV are most probably from cosmic-ray electron bremsstrahlung 
with no associated v-production. The curve marked PCR (|JL ) is an upper limit 
on the v-flux from primordial (or "cosmological") cosmic-ray interactions 

obtained under the assumption that the extragalactic y-ray background 

“6 -2 “1 -l —3 

(8x10 cm s sr above 100 MeV with an E power law spectrum) is from 

interactions of cosmic-rays at high redshifts (Stecker 1969) 0 A detailed 
comparison of the calculated and observed spectra (Stecker 1974, 1978b, 

Montmerle 1977) indicates that the y-ray background is most probably not 
from PCRs so that the curve shown is an upper limit. Note that this upper 
limit is many orders of magnitude below the cosmological neutrino upper limit 
given by Berezinsky and Zatsepin (1977 and references therein). 

The curves in Figure 4 marked H(ATM) and V (ATM) are the horizontal 
and vertical atmospheric background fluxes given by Osborne, Said and 
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Wolf endale (1965) , which indicate the background problems involved in detecting 
cosmic neutrinos. 

The integral fluxes in Figure 5 show galactic "pp" fluxes for the central 

region 330*steA0° (GC) for ]b|£2* and jbjslO 0 , the anticenter region (AC) 

for jb |s2® and the region of the galactic poles (GP). All of these spectra 

were normalized using the SAS-2 y-ray data as presented by Fichtel et al. 

(1978) except for the GP curve which was obtained from Fig. 3 using a mean 

20 2 

galactic path length of hydrogen gas <nL> ■ 3x10 cm . The atmospheric flux 

shown in Fig. 5 is from Zatsepin and Kuzmin (1962). The insert shows the 

galactic latitude distribution of neutrinos calculated for energies above 

1 TeV. TTie curve marked "other galaxies" is obtained by assuming a density 

—76 —3 

of normal galaxies of 3x10 cm (Felten 1977 and private communication) and 

assuming that they have the same y-ray liminosity as our galaxy which is 

assumed to radiate 1/2 of its y-ray flux from TT 0 -decay. Note that since 

the predicted y-ray spectrum from normal galaxies is considerably flatter 

than the observed extragalactic background, so that the limit on the y-ray 

—6 -2 —1 —1 

flux from other galaxies due to tr°-decay is probably ~ 2x10 cm s sr 
above 100 MeV, the v-background from Tt°-decay cannot be significantly higher 
than the curve marked "other galaxies". 

4. Event Rates 

—38 2 

The vN cross section rises linearly with E v and is ~ 0.7x10 E v cm 
at accelerator energies (Barish 1978). The vN cross sections are ~ 1/3-0.4 
of this value (see Table 1), although there is now evidence of an increase 
in this ratio to ~ 0,6 at E v ~ 100 GeV (Barish 1978), possibly due to a 
new flavor b-quark with ~ 5 GeV (Barnett 1976). Above a critical energy 

E v c ■ v^v (1) 
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where My is the ness or the intermediatee vector boson and M is the proton 

w p 

mass, the energy dependence of the cross section levels off to a logarithmic 

one (Bjorken and Paschos 1970). In the unified gauge theory of Weinberg and 

Salem, My * 37,3 esc 8 w where 0^ is a parameter of the theory (Weinberg 1974), 

2 

The best experimentally determined value at present is sin 0 *- 0.3 (Bar is h 

w 

1978) giving My «• 70 GeV „ The resulting yN cross section as a function of 


Ey, is shown in Fig, 6, 

Using the vN cross section from Fig, 6 for My ■ 70 GeV (solid line) 
and a strictly linear cross section (My ■ 00 dashed line) together with 
the calculated cosmic v-f luxes, expected event rates for v^'s falling on 

3 

a 1 km seawater detector were calculated and are given in Fig. 7. The 

dotted line marks the level corresponding to 1 event/yr. The upper 

"pp M curve is for the flux from the galactic central region using the 

production spectrum of Fig. 3 which is derived from the upper v-curve of 

Fig. i. The curve marked is for the galactic "pp" neutrinos averaged 

over 4 tt sr. Thus, on the average, one would expect an event rate of ~1 

yr 1 from galactic v^*s of energy >10 TeV. If ultra high energy cosmic rays 

are universal (not likely) and the vN cross section rises linearly up to 
9 4 

«w 10 GeV (My ^ 5x10 GeV, again not likely) this rate could again be reached 
9 

for Eyj > 10 GeV from yp interactions. It is more likely, however, that 
at energies high enough to see steady diffuse fluxes of cosmic neutrinos 

3 

with a 1 km detector above the atmospheric background, the event rates will 
be very low, 

5. Comparison With Previous Calculations 

Figure 8 shows the predicted fluxes of both and coming from "pp" 
interactions in the galactic anticenter direction as well as the flux of 
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from the yp process* The curves are shown for the present work and 

the results of Silberberg and Shapiro (1977) and Margolis et al. (1978). 

The calculations of Berezinsky and Zatsepin (1977), when adjusted for the 

anticenter region, would be similar to the results of Silberberg and 

Shapiro (1977). The results obtained here, as shown In Fig. 8, give pp 

fluxes which are roughly an order of magnitude below those obtained by 

2 7 

Shapiro and Silberberg for the energy range 10 GeV^Ey^lO GeV and which are, 

on the average, roughly two orders of magnitude below those given by 

Margolis et al. over the energy range 10 GeV^E^e^O GeV for galactic "pp” 

neutrinos. The yp flux calculated here is two orders of magnitude below 

6 8 

that fciven by Margolis et al. for the 10 -10 GeV energy range, (See Appendix II). 

There are three probable reasons why the fluxes given here are lower 
than those given previously by other authors. 

1) Somewhat steeper, but more accurate, primary spectra were used 
in the present calculation. 

2) An erroneous flattening appears in the pp results of Margolis 

2 4 

et al, (1978) in the 10 -10 GeV energy range. (See Appendix I). 

3) Perhaps most importantly, our pp neutrino fluxes were normalized 
to the 100 MeV y-ray observations of the Galaxy, making use of the 
generic relation between charged and neutral pion production (see Fig. 1). 

Finally, it should be noted that existing observations of the extra- 
galactic y-ray background flux and spectrum place severe limits on cosmological 
cosmic-ray pion production (Stecker and Silk 1969, Stecker 1971) and there- 
fore call into question high redshift production burst models such as those 
discussed by Berezinsky and Smirnov (1975) and Berezinsky and Zatsepin (1977). 

In particular, the upper limits given by these authors for the cosmological 
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pp and py v- f luxes are too high. The "pp H upper limit Is that given in 

Figs. 4 and 5 marked PCR(UL) (see Section 3). The py upper limit is not 

related to the X-ray background as proposed by Berezinsky and Smirnov (1975) 

because the X-ray background cannot be related to photopion production as 

originally proposed by Strong et al. (1973) since the shape of the y-ray 

-3 

background spectrum above 10 MeV has the form of «%£ whereas the model of 

-2 

Strong et al would give a spectrum of the form .-JS a Furthermore, the 
Hillas (1967) model on which these cosmological calculations are based, 
implies a y-ray background flux above that observed unless the mean inter- 
galactic gas density at present is <10 ^cm ^ (Stecker 1971, 1975)^, well 
below that implied by recent x-ray observations of galactic clusters, and 
unless galaxy formation occurred at an extremely efficient rate at redshifts 
z^ 15, an assumption which creates problems of compatibility with radio and 
optical observations (Rees 1972). 

The flux of v*s from other normal galaxies is estimated here in Fig. 5 
and discussed in Section 3. Further discussion of possible extragalactic 
v-flur.es is givon by Silberberg and Shapiro (1977). 


^The earlier limit of 10 
given limit of 10 - l°cm - ^ 
Fichtel et al (1978). 


cm was obtained using older data. The presently 
is obtained using the results recently given bv 
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Table 1. 


Ratio 


< V A V 


Vv 


Vv 

/a. 


v e N v e N 

♦valence quarks dominant 

(Eichten et al. 1973a, b). 


yN Cross Section Ratios 


Experimental Value 


1.26 +0.23 
1.32 +0.32 
0.40 +0.12 
0.38 +0.02 


Theoretical Value 
1 

1 

-1/3* 

-1/3* 
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Fig. 1 


Fig. 2 

Fig. 3 

Fig. 4 
Fig. 5 

Fig. 6 

Fig. 
Fig. 


Figure Captions 

. Differential production spectra of neutrinos ahd y-rays from the 
decay of pions produced by interactions of cosmic-rays in our 
galactic neighborhood with interstellar gas having a mean 

3 

hydrogen density of 1 atom per cm . The y-ray curve and the 
upper neutrino curve are calculated for cosmic rays having a 
spectral index of 2.67 between 10 and 3xl0^GeV; The lower neutrino 
curve is for a cosmic ray spectrum with index 2.75. The spread 
in the curve is indicative of the uncertainty in such calculations. 
. Ca 1 nlated differential neutrino production spectra and background 
flu^as from photopion production by ultrahigh energy cosmic rays 
interacting with the 3K universal microwave background radiation. 

. The integral v-production spectrum obtained from the upper neutrino 
curve in Fig. 1. The high energy cosmic ray spectrum used is from 
a recent analysis of Hillas (private communication). 

. Differential background v- fluxes. 

. Integral v~fluxes and the galactic latitude distribution of 
neutrinos with energy above 1 TeV in the direction of the inner 
galaxy. 

Cross section for vN interactions assuming an intermediate 
vector boson mass M w ■ 70 GeV. 

3 

Neutrino event rates for a detector utilizing 1 km of water. 

1. Comparison of present predicted neutrino fluxes with previous 
calculations by other workers. 
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Appendix I 


High Energy Behavior of the Pion-Decay y-ray *nd Neutrino Spectrum 


Let us assume a high energy primary cosmic-ray spectrum of the power- 
lav form 


I (E ) - K E 1 (Al) 

P P P 

interacting with cosmic gas having a density n and producing secondaries 
observed along a llne-of-sight of effective length L. 

We further assume an effective production cross section for pions 
of average energy ^ produced by primaries of energy Ep 

°«nlV 6( VVp b) (A2> 

where ^ and X q are constants. Such a production function would occur 

for a pion multiplicity Q of the form C(E p ) “ The secondary pion 

spectrum is proportional to the Integral over dEp of the product 

cxCeJe )I(E). 

TT P P 

In the high energy limit , y-rays and neutrinos from pion decay would 
then have a power law form with a spectral index 

r- . <r*>- (»+!)- (a3) 

b 

In the case of y-rays, the spectrum would be 

r*3*i 

I(E V ) - 2<hL>K< ^o c ( ~V^ ) E -[ (r+b)- (a+1) ]/b < 

T (r+b)- (a+1) 0 Y 

(Stecker 1971). The corresponding expression for neutrinos is somewhat 
more complicated. 

Again, in the high energy limit, some constant fraction of the 
primary energy E p goes into pions. This fraction is typically of the 


A2 


order of 50Z-60Z, but In any case, It cannot exceed unity. It then 
follows that 

i 

b - 1-a, Oiail, Dibjd. (A5) 

Thus, equation (A3) reduces to 

i' - xEr ■ r + fe (r ' 2 > (A6) 

In order for the total energy in cosir, ic rays to be finite, in 
the high energy limit T must be >2. In fact, observationally this 
is the case. Also, since Q£a£l, it follows that 


f— (r-2)*o, thus r'ir 
l-a 


(A7) 


i.e., the secondary spectrvsn has a spectral index which exceeds or 
equals that of t he primaries. 

In the case a ■ 0, b » 1 corresponding to the dominance of a 

leading pion or scaling, we find r*“r and the spectral index of the 

secondaries equals that of the primaries. It follows from eq. (A7) that 

the spectral index of the secondaries should not be smaller than that 

of the primaries as Margolis, et al. obtain in their scaling model 

2 A 

calculation for the energy range 10 -10 GeV. 

It is well known that if tt* mesons are generated with a spectrum 
g(E n ), the y-ray spectrum from their decay has the form 


g« Y > 


f „ 2 dE n 


f d^g^nl 

Ett 


(A8) 


(e.g. Stecker 1971) in the high energy limit Ey» ra^. Thus if g(E n ) - 


A3 


K B E n ' r n. 8« y ) - K^Y 

*">«* \ ■ V r n- r Y - r„. (A») 

It can also be easily demonstrated that the neutrinos from 
nr -* *iV decay also have the same spectral index as the primaries. In 


this case, the cm energy of the neutrino is 

E v* - = 2 ^ - 29 - 8 Mev - 
Defining 7) * E v */®n» follows that in the limit 


(AID) 


g(E v ) 


2T1 


r ^ 

V21) 


8(En) . JL f* dE _ 8<M 

E" 2 ^ 2.34Ey &n 


(All) 


Again, for a power law pion spectrum of the form g(E n ) * 

— P P 

Vn "> the neutrino spectrum is of the form g(E v ) » ^Ey” v where 


K . (2.3Af r n- 1 ) 


V i\ 




(A12) 


The tt*uv decay kinematics is similar to the Z*‘+Ay decay kinematics 
discussed in Stecker (1971) . The 3-body pr»ew decay is more complicated 
but the result is also r v “ T n (see again, Stecker 1971). 


xnrrmTrn 
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Appendix II 

Approximate yp Neutrino Source Function 

The main source of ultrahlgh energy neutrinos is photomeson produc- 
tion Interactions between ultrahlgh energy cosmic rays and blackbody 
photons of the 2.7K microwave background radiation - i.e., Interactions 
of the type y + p-ei+n. To check on the numerical calculations presented 
in this paper, it is useful to derive an approximate source function for 
the production of neutrinos in photomeson interactions with microwave 
blackbody photons. We a sume all the photons to be at the average energy 
2, 7kTM>.4xlO~^eV so that 

n bb (e) ’ °bb 6(e * e o ) * n^oocm" 3 . (B1) 

The energy of the photon in the cosmic-ray proton rest system is 


e' - (E /M )e rt (l-cose). 

p p O 


(B2) 


There is a large peak in the photomeson production cross-section 

at e'»0.35 M due to the A (12 38) resonance and since f.he cosmic-ray 
P 

spectrum drops off rapidly with increasing energy, most of the pion 
production occurs at this resonance energy. Thus, from Equation (B2) 
we make the approximation 


cr(e' ) a Q 6 


x - 


0.35 




EE 

op 


H Cq6(x 


W- 


where 


(B3) 


X 2 1 - cos0, Eq 
and 

— 2ft 2 
Oq a 2x10 cm . 


~. 35 b. 4.8xl0 ll GeV 
CO 


B2 


Then the source function for neutrino production from THp,v decay may 
be written in the form 


q(E v ) - W2T1 f 


n 


bb 


2 x 

; x dx J dE I(E )o(Ep,x)f (E^jEp) . 
0 0 p p 


We assume all the pions to be produced at the average energy 
<E >- W* 2C-M P . 

" 2 Mp 2 + 2e'Mp 

If we use Equation (B4), Equation (B5) further reduces to 



so that the distribution function is approximated by 


(B5) 


(B6) 


f (E IE ) - 6(E -E /5) . 
tt * p n p 


(B7) 


If we further specify the differential cosmic-ray proton spectrum 

_r 

by a power- law of the form I(Ep) * , Equation (B4) reduces to 

J^ /2T d V E TT 


(E 0 /2) 


^o'bbVo'l J* 


dE E" (r+2) (B8) 


max[E /2,5E V /2T]] 


The solution to Equation (B8) may then be written in the simple form 

. rV 


" (E v> - t Q ^ Ev / El )-(r+i), 


(B9) 


10 , 


where E, s5E /T] =* 2x10 GeV and 
i o • 

Q v “ 2.35xl0~ 23 Kp(E o /2)' r /(m) 


(BIO) 


B3 


For example, if l(Ep) ■ 2.5x10 ^E^ (Linsley 1963) 

Q v = 4.2xl0 _55 cm _3 s" 1 GeV~ 1 . (Bll) 

Equations (B9) and (Bll) may be directly compared with the upper 
curve in figure 2 as a check on the order-of-magnitude of the calcula- 
tions since our numerical results differ significantly from those of 
Margolis et al. (1978). The numerical results shown in Fig. 2 include 
the effect of meson production above the A(1236) resonance. This has 
the effect of producing a flatter spectrum than E v v * at the highest 
energies, one more like E . (Please note that eqs . (8)-(10) of 
Stecker (1973) contain an error in the spectral index but the effect on 
the results of that paper are not significant as indicated by the more 
accurate numerical results of Stecker and Bredekamp (1972)). 


